

protest<-c(0.075, 0.072, 0.072, 0.064, 0.073, 0.069, 0.086, 0.069, 0.065, 0.077, 0.05, 0.066, 0.073, 0.085, 0.075, 0.075)
protest.se<-c(0.028, 0.027, 0.029, 0.025, 0.029, 0.027, 0.024, 0.027, 0.027, 0.027, 0.022, 0.027, 0.029, 0.025, 0.028, 0.027)

var.names<-c("16", "15", "14","13","12","11","10","9","8","7","6","5","4","3","2","1")

y.axis<-length(var.names):1
adjust<-.1


par(mar=c(5,5,2,2),lheight=.5)
plot(protest, y.axis+adjust, type="p", axes=F, xlab="Treatment Effect of Power of Faith Message on Willingness to Protest",ylab="Session Dropped",pch=19,col="black",cex=1.2,xlim=c(-0.15,.15),ylim=c(min(y.axis),max(y.axis+adjust)),main="")
axis(1,pretty(c(-0.15,0.15),5))
axis(2,at=y.axis+adjust,label=var.names,las=1,tick=T,cex.axis=0.75)
box=(bty="1")
segments(protest-qnorm(.95)*protest.se,y.axis+adjust,protest+qnorm(.95)*protest.se,y.axis+adjust,col="black",lwd=1.3)
abline(v=0,lty=2)
